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Abstract 

Three-dimensional source images for mid-rapidity, low transverse momen- 
tum kaon and pion pairs have been extracted from central Au+Au collisions 
data at ^Js^n = 200 GeV by the STAR experiment at RHIC. The pion 
source function displays significant non- Gaussian features implying a finite 
pion emission duration. On the other hand, the kaon source function is es- 
sentially Gaussian, consistent with instantaneous emission from the fireball. 



1 Introduction 

Two-particle interferometry is widely believed to be the most sensitive tech- 
nique available to probe the space-time characteristics of the hot expanding 
fireball created in heavy ion experiments [1]. The first extraction of the 
detailed shape of the 3D source function for pion pairs by the PHENIX ex- 
periment has permitted the decoupling of the spatial and temporal aspects 
of the expanding fireball as well as the determination of a non-vanishing 
pion emission duration using detailed model comparisons [2]. However, a 
consistent picture of the underlying dynamics driving the expansion can be 
achieved only via systematic consideration of a variety of observables. Hence, 
this paper presents the space-time extent of the fireball as obtained from kaon 
interferometry and source imaging from STAR data under similar conditions 
as the pion source extracted by the PHENIX Collaboration. 

2 Experimental Setup and Data Analysis 

The data presented here were taken by the STAR Collaboration during the 
year-2004 and 2007 runs. The colliding beams {y/^NN = 200 GeV) were pro- 
vided by the RHIC accelerator. Charged tracks were detected in the STAR 
Time Projection Chamber (TPC) j3j, covering ±1.1 units of pseudo-rapidity. 
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surrounded by a solenoidal magnet providing a nearly uniform magnetic field 
of 0.5 Tesla along the beam direction. The TPC provides tracking informa- 
tion for track reconstruction and particle identification by means of energy 
loss in the TPC gas. 

The 3D correlation function, C(q), was calculated as the ratio of fore- 
ground to background distributions in relative momentum q for tt^tt^ and 
TT^TT^ pairs as well as K~^K~^ and K^K'^ pairs. Here, q = ^p^) is half 
of the relative momentum between the two particles in the pair CM. sys- 
tem (PCMS) frame. The foreground distribution was obtained using pairs 
of particles from the same event and the background was obtained by pair- 
ing particles from different events. Track merging and splitting effects were 
removed by appropriate cuts in the relevant coordinate space on both the 
foreground and background distributions [4J. There was no significant effect 
on the correlation function due to the momentum resolution of 1% for the 
low transverse momentum under consideration. 

In the Cartesian surface-spherical harmonic decomposition technique [5], 
the 3D correlation function is expressed as 

C(q) - 1 = i?(q) = E E (?X....,(^q) (1) 

where / = 0,1,2,..., = x,y or z, ^ai...a;(^q) Cartesian harmonic 
basis elements (f2q is solid angle in q space) and -Rq^ . q,; (q') are Cartesian 
correlation moments given by 

(^) = J J^q) ^(q) (2) 

The coordinate axes are oriented so that z is parallel to the beam (long) 
direction, x points in the direction of the total momentum of the pair in the 
Locally Co-Moving System (out) and y (side) is perpendicular to x and z. 

The correlation moments, for each order /, can be calculated from the 
measured 3D correlation function using Eq. (|2]). In this analysis, Eq. ([1]) 
is truncated at / = 6 for pions and at / = 4 for kaons and expressed in 
terms of independent moments only. Higher order moments were found to 
be negligible. Up to order 6, there are 10 independent moments: i?^.2, 
Rti, Ryi^ Rt2y2^ ^y6' ^x4y2 and Rl^yi wherc Rl^ is shorthand for 
R^^ etc. These independent moments were extracted as a function of g, by 
fitting the truncated series to the measured 3D correlation function with the 
moments as the parameters of the fit. 

Each independent correlation moment is then imaged using the ID Source 
Imaging code of Brown and Danielewicz [6] to obtain the corresponding 
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source moment for each order /. Bose-Einstein symmetrisation and Coulomb 
interaction (the sources of the observed correlations) are contained in the 
source imaging code. Thereafter, the total source function is constructed by 
combining the source moments for each / as in Eq. ([3]) 

5(r) = E E (^r) (3) 

/ ai...ai 

Alternatively, the 3D source function can also be extracted by directly 
fitting the 3D correlation function with an assumed functional form for S{t). 
This corresponds to a simultaneous fit of the independent moments with the 
assumed functional form via the 3D Koonin-Pratt equation Eq. (jl]) 

C(q) - 1 = Riq) = J rfri^(q, r)S(r) (4) 

3 Results 

Figure [T] (a) compares the / = moment (solid circle) with the ID correla- 
tion function (open circle), for mid-rapidity {\y\ < 0.35, where y is particle 
rapidity), low tranverse momentum kr (0.2 < kx < 0.36 GeV/c) 7r+7r+ and 
TT^TT^pairs from 200 AGeV central (20%) Au+Au collisions, as a function of 
the relative momentum q in the PCMS. They are in very close agreement, as 
is expected in the absence of significant angular acceptance issues. 

The model-independent imaging technique of Brown and Danielewicz [6] 
is used to image each relevant correlation moment to obtain the correspond- 
ing source moment. The latter are then combined according to Eq. [3] to yield 
the 3D source function. 

Figured] (b), (c) and (d) compares the magnitude of the correlation mo- 
ments (symbols), for increasing orders /, in the x, y and z directions respec- 
tively. The moments decrease in magnitude with increasing order and are 
virtually for Z = 6. This justifies truncating the series Eq. ([1]) at / = 6. 

The restored correlation moments, obtained from imaging the actual cor- 
relation moments, are shown by lines in Figure [T] (b)-(d) for I = 0,2,4,6. 
They are in good agreement with the actual correlation moments, implying 
that the imaging process was well under control. 

Figure [2] (a)-(c) show profiles of the extracted 3D source function in the 
X, y and z directions obtained by summing the source moments up to / = 6. 
The source image (stars) is in good agreement with that reported by the 
PHENIX Collaboration (squares) in all 3 directions [2]. Panels (d)-(f) show 
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Figure 1: Pion correlation moments (open symbols) for multipolarity orders (a) 
1=0 (b) 1=2 (c) 1=4 and (d) 1=6. Restored moments are shown as curves. Panel 
(a) also shows a comparison between R^{q) and R{q) = C{q) — 1. 



that the corresponding correlation function profiles from STAR and PHENIX 
are also in good agreement. One can also note that the broader source profile 
in the x direction is associated with the narrower correlation profile. 

Figured (a)-(f) (left) shows the correlation moments (open circles) ob- 
tained from the cartesian surface-spherical harmonic decomposition of the 
3D correlation function for mid-rapidity {\y\ < 0.5), low kx (0.2 < < 
0.36 GeV/c) K+K+ and K-R-paiis from 200 AGeV central (20%) colli- 
sions. Panel (a) shows that the / = moment (open circle) is in good 
agreement with the ID correlation function (solid circle). The / = 2 mo- 
ments in panels (b)-(c) are relatively small and the I = 4 moments in panels 
(d)-(f) are statistically consistent with 0. The curve represents an ellipsoid 
fit to the six independent moments. The /ndf value of 1.7 indicates that 
the kaon source function is reasonably described by an ellipsoidal shape. 

The dotted curve in figure[3](a)-(c) (right) shows the profile of the ellipsoid 
fit in the x, y and z directions. The band around the curve indicates the 
statistical uncertainty on the source fit. The extracted kaon source function 
is fairly well reproduced by the Therminator model [7] calculation in blast- 
wave mode (triangles) with a transverse size of pmax = 7.5 fm and lifetime 
r = 7.5 fm/c of the emitting fireball and kaon emission duration Ar = 0. 
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Figure 2: Comparison of pion source function profiles S{rx), S{ry) and S{rz) (pan- 
els (a)-(c)) and their associated correlation profiles C{qx), C{qy) and 0(^2) (panels 
(d)-(f)) in the PCMS from the STAR (stars) and PHENIX (squares) experiments. 
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Figure 3: Panels (a)-(f) left: Correlation moments i?'(g) for I 
also shows a comparison between R^{q) and ID R{q). The solid lines indicate the 
ellipsoid fit. Panels (a)-(c) right: Kaon source function profiles from STAR. Solid 
triangles depict Therminator model calculation with indicated parameter values. 
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4 Discussions 



The difference in the extracted 3D source function between pions and kaons 
coupled with the determination of similar spatio-temporal characteristics of 
the expanding fireball from both pion and kaon correlation functions suggests 
that pion femtoscopy is significantly influenced by resonance decays and emis- 
sion duration effects compared to that for kaon. Indeed, a non- vanishing pion 
emission duration of about 2 fm/c was reported by the PHENIX Collabo- 
ration |2] whereas the new STAR result points to kaon emission as being 
essentially instantaneous within the framework of the Therminator model . 
Furthermore, the kinetic freeze-out process for pions is different from that 
for kaons as evidenced by the different values for the slope parameter, a, 
of the freeze-out hypersurface. The latter specifies how the freeze-out time 
varies with the transverse radius of the expanding fireball. While a negative 
value for a was necessary for pion emission, implying emission of pions from 
outside-in, a zero value for a is sufficient here to explain kaon emission. 
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